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ABSTRACT 

The major chemical systems of the Earth’s atmosphere are tied to absorp- 
tion of radiation from the sun, and thus have a dawn and dusk dependence. Also, 
total solar eclipses may cause sharp local variations. 

A study is made of the behavior of chemically generated waves in 
a simplified atmosphere. The atmosphere is assumed unbounded, isothermal, 
one-space-dimensional and initially quiescent. At an initial time a dissociation 
reaction commences and drives the subsequent wave motion. The fraction of 
reactant in the atmosphere, X, is assumed to be small. The system of govern- 
ing equations is then expanded in terms of the small parameter, X, and an in- 
tegral solution asymptotic to X -• 0 is obtained. A series solution in terms of 
Young functions is also obtained. These are a subspecies of Lommel functions. 


♦Now at the Division of Applied Mathematics, U. S. Naval Ordnance Laboratory. 
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A LINEARIZED APPROACH TO WAVES 
IN A DILUTE, CHEMICALLY REACTING ATMOSPHERE 

PART I 

ANALYTICAL TREATMENT 


INTRODUCTION 

Waves in a stratified fluid under the influence of gravity appear to have 
been initially discussed by Burnside (1889) and Love (1891). Both authors 
treated an incompressible fluid. Gortler (1943) used schlierin photography to 
show experimentally that disturbances in an incompressible stratified medium 
under the influence of gravity propagate along characteristic rays. Lamb (1908) 
treated a compressible, adiabatic, ideal gas whose density is stratified by 
gravity. Since the original contributions by Burnside, Love and Lamb there 
have been many theoretical papers on various types of gravity waves. Reason- 
ably up-to-date treatments of the subject may be found in Eckart (1960) and Yih 
(1965). However, gravity wave theory is currently undergoing an active phase 
of development and contributions are being added to the literature at a rapid 
rate. 

Recent experimental data from the earth’s atmosphere and oceans have 
greatly stimulated interest in gravity-acoustic waves. Gossard (1962) has ob- 
served gravity waves in the troposphere. Hines (1960) has shown that in- 
ternal atmospheric gravity waves may account for many of the phenomena 
observed in the lower ionosphere, and gravity waves have been observed in the 
thermosphere by Newton et. al., (1969) and Harris et. al., (1969). 

It seems to be generally agreed that gravity waves are generated in the 
troposphere, then propagate upwards. Thus Gossard (1962) has observed gravity 
waves near the earth’s surface and Eber stein (1970) has illustrated the develop- 
ment *of gravity waves between 30 km and 120 km. 

The simpler forms of gravity wave theory treat propagation through an 
inviscid, non-conducting gas, while more sophisticated theories include effects 
of heat transfer, viscosity, and high altitude phenomena such as ion drag. How- 
ever, the effect of chemical reactions does not seem to have been adequately 
considered to date. 

Between the troposphere and the thermosphere there are several regions 
where important chemical reactions take place. 


1 


First, there is the ozonosphere between approximately 15 km and 45 km, 
with peak ozone concentration at about 35 km. Ozone mole fraction is in the 
order of parts per million (Mitra, 1952). 

Between 60 km and 96 km there is Nitric Oxide (Pearce, 1969). The mixing 
ratio for Nitric Oxide is also in ppm. 

Between 90 km and 120 km molecular oxygen dissociates into atomic oxygen. 
The dissociation ratio, a being 3 x 10" 8 at 90 km and 0.998 at 120 km (Mitra, 1952). 
The oxygen is no longer a truly dilute reactant since its mixing ratio is 20%. 

The major chemical systems are tied to the absorption of ultraviolet radia- 
tion from the sun, and are thus have a dawn and dusk dependence. 

In what follows, an initial study is made of atmosphere. Simplifying 
assumptions are made to make the mathematics more tractable. 


2. Governing Equations 

The continuity and momentum equations (cf. Shere and Bowhill, 1969) are 
not affected by the chemical reaction provided that the state variables are given 
the extended definition: 


P . PO, T, a) 


where P, p 9 T and a denote pressure, density, temperature and degree of advance- 
ment of the reaction, respectively. These equations are in tensor notation: 


and 


Sr 

DT 


+ P 


du i 

dx. 


0 


( 2 . 1 ) 


a 1X1 i _<9P a <7j j 

Dt ~~ P% i ~ dx A + dx. ( 1 


( 2 . 2 ) 
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where 


du i du. 2 3u k ‘ 

a H a ^[ 57 : + 377 •T s ij75^_ 


We denote time, space variable (x 2 vertical), velocity component, gravitation 
vector component, viscous stress tensor, dynamic viscosity and the Kronecker 
delta by t, x., u. # g i# a. ., fi and 3^ respectively. We also take g* “ “gS i2 and 
we define D/Dt to be the total or Euler ian derivative. 

In order to develop the energy equation, it is necessary to consider the 
thermodynamics of the system. The derivation of the equation of state parallels 
the work of Iiepmann and Roshko (1957, p. 29) and the derivation of the reaction 
rate equation parallels the work of Eberstein (1966). The following discussion 
applies to the reaction 


J + AB 71 A + B + J 


where J is an inert element. We assume that each component of the mixture is 
a perfect gas. The mixture, however, is not a perfect gas. The partial pres- 
sure of the i th constituent is 


P. = n>. P R. T 


(2.3) 


where R. = R/W i with R the universal gas constant, W L the molecular weight, 
and m. the mass fraction of i th- component. Letting a be the fraction of AB dis- 
sociated, n i the number of moles of ith component, and n 0 the total number of 
moles when a = 0 , 


AB 

✓—s 

a 

I 

H 

O 

X 

0 

c 

II 


n A 

n o *o a 

(2.4) 

n B 

~ n o ^0 a 


n i 

= (1 - X 0 )n 0 
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where X 0 is the mole fraction of AB when a = 0. Summing the partial pressures 
yields 


P = R(a)pT 


(2.5) 


where 


R(a) = R(0) ( 1 + X Q a) 


( 2 . 6 ) 


and 


R(0) = R/[Wj + X 0 (W AB - Wj,)] 


The reaction equation is as follows: 


DN AB /Dt 



(2.7) 


where 



and k is a rate constant. In determining (2.7) it has been implicity assumed that 
all molecules present may act as an inert element with equal collision effectiveness. 
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Substituting (2.4) into (2.7) yields 


Da 


Dt = [pk,(T)/M] 


Jl - a 


l 


[K(T)p/M] 


(jy_ 

\ll + X 0 a] 


(2.8) 


where 


t 

V 

*/# 


k„(T) 

k(T) . crhry 


r* 

v 

« 

V 

. 

:$■ 


and M is the mean molecular weight 


M = + M*ab -«j> 


For chemical equilibrium, the rate equation reduces to 


-a. ± (-**-) 

— K p \l+axJ 


where use was made of the equality 


K = 


Kp 

K c " RT 


lh atmospheric problems the reactant fraction is generally small, and it 
will be assumed that X 0 < <1. It follows that aX 0 « 1. The rate equation may 
thus be simplified to give: 


Da 


5F = [pk p (T)/M] {l - a - [K(T)p/M] X 0 a 2 } 


(2.9) 
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The energy equation is now developed in the usual manner (cf. Liepmann 
and Roshko, 1957, p. 185+). Summing the rate of change of the internal, kinetic 
and potential energy yields 


Jjt (-° e + T*»i u i + P%*i) dV 


J + -j/OUiUi + u j n j dA 


- J Pn t Uj dA + Ju.cr..dA+ Jk ^ n. dA < 2 - 10 ) 


y are n. is the ith component of the unit normal of the closed surface <£, V is 
the volume of and A the area; e is the internal energy per unit mass and k is 
the coefficient of thermal conductivity. 

By applying Gauss’ theorem and substituting the continuity and momentum 
equations into (2.10), we obtain 


De du i du. L d /, JI\ 

^"Dt = " p a77 + + [ k 9xJ (2.ii) 


We now use the first law of thermodynamics, h - e + p/p, and the continuity 
equation to obtain 


Dh _ DP 
^Bt Dt 


du 


+ cr. 


» i dx . 



( 2 . 12 ) 


In a reacting gas mixture the enthalpy depends on the temperature and the 
degree of reaction, i.e. 


h = h (T, a) 


The rate of change of enthalpy per unit mass is then given by: 


Dh 

T5F = 


DT 

r 

Dt 


+ X a B 


Da 

Dt 


(2.13) 


where B is the enthalpy of reaction and C p is heat capacity. Both the enthalpy of 
reaction and the heat capacity have a weak temperature dependence, and an even 
weaker pressure dependence. For simplicity both C p and B will be assumed 
constant throughout this paper. 

Eliminating h between (2.12) and (2.13) one obtains the enthalpy equation in 
the form: 


pCp-i 


DT 

Dt 


X„B 


Da 

Dt 


DP fill , JL /. _iX\ 
Dt + CT ik <9x k + «?x. Sx J 


<2.14) 


To summarize, we have found the system 




0 


( 2 . 1 ) 


D^ 
P Dt 



tin. 

dx, 


( 2 . 2 ) 


P = pR(a)T 


(2.5) 
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« ■ ■ ' V-isw ■ ; « ~iTi . vi ,i i i ;ii, v -- i.j y.v > — , „n . 


■ ^ a a a ar 


Da 

Dt 



DT Da DP 

(G . p Dt + X o P B Dt = Dt 


K(T)X„ a* 
-a - p h 

m§ 


^ i j_ ( 3T \ 
+ a H 5Xj + axj \ k <?xJ 


(2.9) 


( 2 . 11 ) 


The dependent variables are P , p, T, u i , and a; the independent variables 
are t, x .. k F , K, c p and B are known parameters of the system, X 0 is a small 
(known) parameter. 

For the remainder of this paper w r e assume that the atmosphere is inviscid 
and nonconducting. Hence, the governing equations are 


Dp 

Dt 


+ 


du. 



0 


(2.15) 


Du. 

i 

P Dt 


<9P 

dx k 


(2.16) 


P - pRT 


(2.17) 



(2.18) 


P C p Dt + X o B P Dt = Dt (2.19) 

3, Nondimensionalization of the Governing Equations 

Since mathematical operations are performed on pure numbers it is proper 
to non-dimensionalize the governing equations. 
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Some of the non-dimensionalization is almost trivial, as shown immediately 
below: 


R' 

= R/R(0) 

V 

= C p /R(0) 

p' 

= P/P , 

M' 

= M/M = 1 

k # 
k F 


K' 

= Kp./M^ 


where R(0) is the gas constant for a « 0. 
p* is a reference density. 

Other non-dimensionalizations have a more direct physical or kinematic 
meaning. Thus, distance is defined in terms of scale height and time is related 
to the Brunt- Vaisala frequency, i. e. 

*1* = *l/H 

t ' = ta> B 

Further non-dimensionalizations are: 


T' = R(0) TA> B 2 H 2 

B' = B/oj b 2 H 2 
g' = g/Hai B 2 

P' = ai B 2 H 2 
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It is seen that in the primed system we have 


P' 


p' R' T' 


In general, the non-dimensional system of governing equations looks very 
similar to the dimensional system. 

The quantity <^ B 2 H 2 is proportional to the energy contained in a wave 
oscillating with the Brunt frequency, and having amplitude equal to a pressure 
scale height. The quantity Ho> B 2 is a measure of the acceleration experienced 
in such a wave. Since such acceleration must be less than that due to gravity 
we would expect 


g' > 1 

One may estimate the magnitudes of the non-dimensional parameters. 
Taking a scale height of 6 km and a Brunt period of 5 minutes, one obtains 


(<^ B H) 2 ~ 0. 1 calories/gram 


and 


CO 


B 


H ~ 5 cm/sec 2 , giving 


g 

T # 


~ 200 
~ 200 


Further, taking B = 10 kcal/gmole one obtains B'~10 s . 
4. Asymptotic Development 


The non-dimensionalized system of equations will be expanded about the 
parameter X 0 , i.e., we expand each dependent variable in a possible series of 
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the form: 

00 

f(t,x, z) = ^ f (N) (t, x, z) X/ < 4>1 > 

N *0 

where higher order terms may be dropped as X 0 -♦ 0. 

Assume that the atmosphere is initially stratified, quiescent, and bounded 
below by the earth. 

Hence we have defined an initial value problem with initial conditions: 

T(0, x, z) = T 0 = T, 
u.(0, x, z) = 0 (i = 1, 2) 

(4.2) 

P( 0, x, z) = e" * 
a(0, x, z) = 0 

and boundary condition: 

u 2 (t, x, 0) = 0 

Using the equation of state to eliminate P from (2.18) - (2.19) and substituting 
(4.1) into the resultant system yields the zero-order solution 

T<°> (t, x, z) = T 0 = T. 

u. (0) (t, X, z) = 0 

(4.3) 

P (0) (t , x, z) = e 2 

a<°> (t, x, z) = 1 “ exp [-k F /3< 0) t] 
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Defining TO) = T<°> T (l) , P<‘>= P<°> P (1 > anda<*> =a«»a (l) . we get the 
first order system 

au.,,. , T 3|0 O> . T aT( = 0 

"57' ) + T o dx o dx (4,4) 


3u 2 <‘> 


dP (l) 

+ T 0 dz ‘(0 ) ‘(i) 


dT 

~ T/ ft v T f i ) + T 0 dz 


i±L = T a (0) (4.5) 


^n + !^ + ^i- u <n = o. 

— - + “ + Jz u 2 


(4.6) 


at dx 


^ a ( 1 ) 

at 


i ^(D 

+ y-=“r at 



= -B 


aa<°) 

at 


(4.7) 


^ a <>) 

at 


a 


( 0 ) 


aa(°> 

+ a <i) at 


r (l) (p<°>, a<°> ) 


(4.8) 


where 


8 = (t 0 - *) 


The above forms a linear system of four partial differential equations in 

four unknown functions. The system is inhomogeneous. The solutio " ( ■ ' } 

for a (1) is needed only for the calculation of second order terms. In this paper 

we determine only the first order terms. 

Dropping the (1) subscript or superscript and eliminating we obtain the 


equation 


Liirii - 

(j3t 2 l_at 2 


c2 il-c 2 4 + c 2 a 
a* 2 az 2 


a “I y~ i 
■5¥j ~~yi~ 


C 4 = f ( a 

ax 2 


(4.9) 
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where 


f(a<°>) * 


-(7-1) 





- (7 - 1) B 


a 3 a(°> 
dt 2 5 z 


d 2 ® ^ 


and 


C 2 ^ y T 0 

Using the results of Garding (1950), A, Lax (1956), and Courant and A. Lax 
(1955) it can be shown that (4.9) is a well-posed hyperbolic equation. 

After (4.9) is solved the other dependent variables are determined by the 
equation 



+ T, 


±] 

dzj 


A. All 

7“i at 2 


dT 


+ T 0 T ' + T o a<0> < 4>11 > 



d 2 P 
d t < 9 z 


A. 32 T 

+ y - 1 3 t dz 


A il 

y- l dt 



( 4 . 12 ) 


and 


_ 3 P 

u 2 " TT 


1 < 9 T „ da<-°'> 

7TTTT ” B Tt 


Now define 

0 (t, x, z) = e ~ t/2 T (t, x, z) 


( 4 . 13 ) 


(4.14) 
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* 


The above definition allows the temperature perturbation to grow exponentially 
with altitude without a corresponding growth in 6 . 

Substituting the definition of 0 into (4.9) yields: 


d 2 
dt 2 



- C 2 V 2 



c 2 (ri) d 2 e 

7 2 <9x 2 


e “*/2 f(a) 


(4.15 


f (a<°>) was defined in Equation (4.10). 

The dispersion equation of (4.15) is: 

k* ■ c * ( X x 2 + k z + t) k t 2 + [(o' ” 1 ) c 4 /y] K 1 = 0 < 4 * 16 > 


where k i is the wave number of i. 

For a discussion of dispersion relations and how these are obtained for 
partial differential equations reference is made to Courant and Hilbert, especially 
p. 588. 

5. Solution of the First-Order System 


t 


In this section we first demonstrate that the governing temperature Equa- 
tion (4.16) is reducible to an inhomogeneous telegraph equation with homogeneous 
initial conditions. For a nonreacting atmosphere which is initially perturbed, 
Lamb [190 9] ; obtained a homogeneous telegraph equation with inhomogeneous 
initial conditions. 


We take for initial conditions of 9 (t , x , z ) , 



0(0, x, z) " 0 t (O, x, z) 

= 0 

(5.1a) 

0 tt (O, x, z) = -(y- 1) 13 k F a 

exp (-5z/2) 

(5.1b) 

6? ttt (0, X, z) = (7-1) [® kp p < °> 2 

+ T 0 ]exp (-3z/2) 

(5.1c) 


14 



Conditions (5.1b) and (5.1c) are determined by the requirement that 0(t, x, z) 
remain bounded as t - 00 . The negativeness of 8 tt implies that the atmosphere 
initially cools. This is in correspondence with the endothermic reaction chosen. 

Since the inhomogeneous parts of (4.15) - (5.1) do not depend upon the hori- 
zontal space variable, neither does the solution; this is obvious from a considera- 
tion of the Fourier representation of <9(t, x, z) . Thus (4.15) reduces to 


d 2 

dt 2 


M[0] 


e z/2 f (a) 


where M C • ] is the telegraph operator 


(5.2) 


MH = <9 2 /d t 2 -C 2 d 2 ,/d z 2 + C 2 /4 . 


(5.3) 


It is worth noting that the initial conditions and the alpha dependence chosen 
are such as to imply the relevance of only one space dimension. Other initial or 
boundary conditions would require that the second space dimension be kept. In- 
tegrating (5.2) with respect to time twice, 

M[0] = w(t,*) = k F (yl) [(1-T 0 t 2 ) Bk f e 1 

+ (213 +1) (C 2 /y)t]e* 3z/2 " k F e ’* < 5 - 4 ) 


subject to initial conditions (5.1a). 

The inhomogeneous part of (5.4) is the effect of the chemical reaction. For 
a realistic approximation to the atmosphere, the inhomogeneous part of (5.4) must 
be multiplied by a suitable weighting function, since the chemical contribution only 
extends over a finite altitude regime. We also note that at each fixed altitude the 
solution must tend asymptotically to zero in the Poincare sense as t -* 00 • 

The assumptions used in this paper are also suitable for a multilayer ap- 
proximation of the atmosphere. In this case we would consider either an initial- 
boundary or pure boundary value problem for (4.15). The x-dependence may not 
generally be eliminated for this case. 
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The solution of (5.4) - (5.1a) is well-known [cf. Courant and Hilbert, II, 
pp. 695, 202] and is given by 


0(t, z) 



W (t, z; t) dr 


(5.5) 


where 


W(t, z; r) 



Q(z , 7 )) drj (5.6) 


and 


Q< z, 17 ) a [w(t, z+ T 7 )-w(r, z - 77 )] /C 2 (5.7) 


Although (5.5) - (5.7) provides an exact solution of (5.4) - (5.1a), this solution 
is complicated and does not easily yield qualitative information. 

The first order terms are from (4.12) - (4.13): 


T(t, z) = e* /2 6( t, z) 


u 2 (t, z) 



where z § is at the earth surface and 



To obtain qualitative information we consider the Laplace transform of (5,4): 


■W„(*, <0 + (<r 2 + l/4)W = (r-1) (2B+l)(k r /7C) 


»• 3 * / 2 


[o- + (k f /C)e‘ , l* 


B(k F /C) 2 (y-l)e 


[<r+ (k r /C)e'*] 


-5t/2 


B(k F /C) s [(r-l)/7]e’ 5,/J 
[or + (k p /C) e" *J 3 


(5.8) 


where W(z, cr) = £|i9(z , r)\ and r = tc . 

Equation (5.8) is most readily solved by superposition. We obtain 


W(z 


0 f (~l) N ( k F^ c ) N ,1 

, CT ) = 2- B(kp/C) (r_1) V^[ CT 2 -(N + 3)(N + 2)]j 

N - 0 J 


> g “(5/2 +N>* 


„ f (n + l) (“k F /C) n 1 

L <r ' 1) < IB *»<vw>L... iTj 

n = 0 v 


L e -(3/2 + n ) * 


“ B(k r /C) a (r- 1) f (n + l) (n + 2>(-f r /0r)»^ 

■ ■ ^ «' 
n* 0 


27 |«r n + 3 ^. a 2 + 3) (n + 2)]J 


v e “ (5/2 \+ n )* 


(5.9) 
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Using Bateman's formula (10) [Erdelyi, et. al. I, pg. 238, 1954] we take the 
inverse Laplace transform: 


0(t, z) 


V" 1 /( ? - l)13(-k F t) n * 2 / n-^3 n + 4 (nt 2 )(n +3)C 2 t 2 \ 

/ [ r (n + 3) I F 2 \ ; 2 » 2*» ^ / 


n = 0 


(y ~ 1) ("*^f *) n * 4 / n + 5 n + 6 (n + 2)(n+3 )C 2 1 2 \ jw - 

+ 2yk_ a T (n + 5) ‘ Fj V 1 ' 2-2= 4 /J e 


(S/2 + n) j 


& 



(y-l)(2l3 + l)C a (-k F t) n + 3 
rk F P(n + 4) 


P 
1 2 



n + 4 
2 ’ 


n + 5 (n41)(nf2)cV , (3/2tn ^ 

2 ’ 4 / e 

(5.10) 


We identify the generalized hypergeometric functions appearing in (5.10) as 
Young's functions [Young, 1912], a special case of Lommel functions. 

It had been noted earlier that £(t, z) must be asymptotic to zero as 
t oo . This behavior in time can be obtained by applying Luke's formulae 
5.11.2 (T) and 5.11.2 (9) [i, 1969, pg. 199] . The series (5.10) is rapidly covergent. 

The results of some calculations on the chemical effect are given in Part II. 
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